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IV. On some Definite Integrals, and a New Method of Reducing a Function o y 

Spherical Co-ordinates to a Series of Spherical Harmonics, 

By Arthur Schuster, FM.S. 



Bocewed Mav 30, — -Read June 5, 1902. 



§ 1. Introductory. 

The following investigation deals with some definite integrals which are useful when 
it is desired to express a function of two angular variables by means of a series of 
spherical surface harmonics. An important theorem concerning these integrals leads 
to a method which considerably reduces the arithmetical labour involved in the 
reductions, and secures in practice the advantage of obtaining the numerical values 
of the coefficients of lower degrees independently of those of higher degrees. 
The zonal harmonic of degree n is denoted by P w and defined as usual by 

1 rl n 

P -= (a 2 — lY 

u " 2* 1-2... ndfS l ^ ; * 

The tesseral harmbnic of degree a and type a is denoted by 

T;I = QZ (gl cos <x<i> + K sin cr<£), 

d°-V n 



where QZ = sin 47 6 

'JjJL 



11 ii-v 



In these equations 6 represents an angle jneasured on a sphere from a point as pole 
and fi = cos 0. The longitude measured from some standard meridian is denoted 
by (j>. 

The name given by Heine, and translated by Todhuxter as i; Associated Functions 
of the First Kind/' is too cumbersome for use, and I propose to call these functions 
" Tesseral Functions." The tesseral function is converted into a tesseral harmonic by 
the factor cos or/> or sin cr</>. The name may not perhaps appear to be appropriate, 
because it is only the factor which o'ives the function its "tesseral" character, but it 
is short and suir.^ests at once the function it denotes. 

The present investigation deals in groat part with the definite integrals, taken over 
the surface of a sphere of unit radius, of the product of two tesseral functions and of 
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the product of a tessera! function and cos p0 or sin p9 9 the tesseral functions being all 
referred to the same axis. I shall denote, as usual, the factorial product of the 
numbers up to n by n ! , but have found it necessary to introduce a separate notation 
for the products of successive even or successive odd numbers. I consequently define 

n \\ = n. (n — 2) ! ! . 
Starting with 1 ! t = 1, 2 ! ! = 



o 



it follows that, for positive values of 



i v 






where the last factor is either 1 or 2, according as n is odd or even. 

We may extend the definition to negative values of the argument, for the successive 
substitution of n = 1, n= — 1, n = — 3, into the first equation, gives 

(_1)!! = 1, (-3)!! = -I, (-5)!! = i 

and generally if n Is negative and odd, 



„!! = (_ ifr 1 



1 






For n ™ 2, the original equation gives ! ! = 1, 
n = 0, (-2)!!= cc, 

and similarly for all negative and even values of w, n I ! Is infinite. The ratio of two 
of these factorials of negative numbers is, however, finite, for it is easily shewn that 

if m and n be two negative numbers, whether even or odd, 

One of the advantages of a separate notation for what may be called the " alternate " 
or " double" factorial, is due to the fact that It often saves the inconvenience of 
different expressions for odd and even numbers. 



| 2. Formulas of Trans formation. 

It is convenient to collect together some equations which will often be required. 
Most of these equations will be found already in previous writings, such as Heine's 
Treatise or Adams' " Researches in. Terrestrial Magnetism." 

As regards zonal harmonics, it is only necessary to quote the well-known relations ; 
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(2n+ I) /*?„ = («+ I)P M+1 + ^P^x (1), 

(sn + i)p.=^r J - v 1 «• 

From these equations we derive the following : 

(2H. + i) / iQ; = (-»i-o- + i)Qj +1 + (n+or)Q;_ 1 (A), 

(2«+l)sm^Q^ = Q^;-Q^i . (B), 

= (n + a) (n + <r - 1) Q^lJ - (n - o- + 2) (w - <r + 1) (&;} . (C), 

gf n 9 f=(w + o-)(n + .o--i)Qi;:I + Q;ii (D), 

= Q?:i + (w-<7 + 2)(n-o-+l)Q^l (E), 

Q;-Q;-2 = (»+o--2)(« + o— 3) Q;:i-(»- o- + 2)0 4 -<r+i)Qr a •■ (F), 

(2ii 4- 1) -■■■- sin'fl -- 

= (n+ o-) (n + o-—p+ 1 ) sm'--0 rf ' I> ";-- 1 - (« - 0-4- 1) (w- 0-4-/0) sin" -°-O d 'fw . (G), 

2(x y- sin" -■ p = (2cr - p) sin' 5 ^f - p (n + cr) (« - a + 1 ) sin'" 2 —---"- . (H). 

As special cases of (G) and (H) we may put p equal successively to cr, or + 1, and 
cr + 2. The following equations are thus obtained : 

(2» 4- I ) sin- 6- j Ql = (/i 4- 1) (n + cr) <&_, - n (n 4 1 - cr) Qj +1 . (G L ), 

7 

2 sin (9 ^ Q: = Q: " - (w 4- cr) (n - a- + 1) Qj" 1 .... (HJ, 

(2« 4 1) sin ^ sin Ql = n (n + a) Qj_, - (n 4- 1) (« 4 1 - cr) Q* +1 . (G 3 ), 

2a- ■}- sin 6 Ql = (cr - 1 ) QT ' - (cr + l) (,i 4- cr) (-« - <r 4 1 ) QT ' • (H 2 ), 

(2u + I) -~ sin- 6? Q; = (/i — 1) (u + a) QU - (n + 2) (n - cr + l) Q; +1 (G 8 ), 

2cr — sin 2 6 Q; = (cr - 2) sin Q;; +1 - (cr 4- 2) (n + a) (n - cr + 1) sin Qr 1 (H 3 ). 

The formula (A) is well known and may be obtained by cr differentiations of (1), 
substituting in the result an equation derived from cr — 1 differentiations of (2). 
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Equation (B) is the result of cr differentiations of (2), and multiplication by 



sin <r+1 9. 



C) may be proved by combining 



d^ K /X ) dfi 



f(T + 1 



(I-/* 2 ) 



70" T 11 J 



rf//, 



/ 3 ~ T > 



<x+l 



a/A 



1 



rZ/x 0- 



with the fundamental equation for zonal harmonics. The latter leads directly to 



dl\ 



cr 

/l 2\ 



cl 



7 rr — 1 

' a ^n(n + I). P 



il 



(JL 



from which we derive (by equating the two expressions on the right-hand sides) : 



d' J+1 l y 



rrp, 



« 



cr 



7 cr 
(tfl 



(u -f cr)(n 



(X 



■Hi)',..- " 



cr 



or, after multiplication by (1 — ju~) 2 , 



sin 0Q£ +1 = 2/xorQ,; — (v/ + cr) (n — cr + 1) sin 0Q*" 1 . 



if /xQ* be now substituted from (A) and. sin^Qr -1 from (B), the equation (C) is 
obtained. 

If cr — I be written for <x in (B) and cr + 1 for cr in ((J) we may combine the two 
equations, so as to give (D) and (E). (F) is an important relation obtained from (B) 
and (0). The formulae (G) and (H) are easily derived by direct differentiation and a 
few simple transformations. 



§ 



+i 



o 



O a dfi. 



%,) -* 



If the equation (BL) is integrated with respect to fx between the limits — 1 and 
+ I, the left-hand side vanishes at both limits: hence, after changing from a + 1 
to o\ 



r + \ 



Ql dfx 



<i 



r-M 



a 



(n + cr — I) ("it — cr + 2) j Q*" 2 dji. 



and, by applying the same process to the right-hand side, 



<T 



a 



4: 



(oi + cr — l) (7^ + ^ -"" 3) ( n 



c + l. 



cr 



}.,, o 



) (-ft — cr -f- 4) Q* 4 dji 



Repetition of the same proceeding will ultimately lead for even values of cr to 



+ 1 



-1 



Qn dp 



cr (n + cr ~~ 1) ! ! (w —2)!! 

2 (^™~a-)!! (7i'+l)!! 



;>\ t » - + l 



Qh ^- 
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But 






(1 — /x 2 ) : 'y~~ du 
4 ' J CIJJ," ' 



M 

AJ 



4- 1 dX> 



2 






P„ j dfi 



2 



/*P. 



+ 1 



where the special case n = 0, for which the integral vanishes, is excluded, 

xxeiivv 

■+l 

Q^ dfi = 4, if >?/ be even ; = 0, if n he odd ; 



W **.-*# J 



and finally 



+i 



"Q; fy 



(n + <r-l)!! 0/4-2)!! .« , , 

2<x ; -;■,, --, ;,;, , n a 1 and 7i be even 

(n — a) ! ! (?i. -f 1) i ! 

0, if a* be even and n odd. 



For odd values of <x, it is more convenient to use formula (G 3 ), which leads to : 






■l 



(?& — cr) (w + 1) 



By repetition of the same process, we get for even values of n ultimately 



+ i 



^■cJjil C''M/ 



n + or — 1) ! ! ('//, 



9"\ ? ? 



cr '. 



•+i 



0-<r) !!(?& + 1)!! (<r-3)!! (2<r - 2) ! ! .Li 



Q*-i ^ 



As the expression under the integral sign of the right-hand side vanishes, the 
value of the integral on the left-hand side is zero. 
When n is odd we are ultimately led to 



+ i 



■/£ --fa* — 1 



t » 



7i 



ON t * 
I) l I 



cr 4- 1 



» * 



O - a-) ! ! (m + 1) ! ! (a- - 2) ! ! (2<r - 1) ! ! 



4 1 

V^g- Cf Lilt 



But 



+ 1 
Qj dp ^ (2or — 1) ! ! 



+ i 



sin 



rr 



dfi = ( 



m ir 



2a 



1) ! ! | sin 

o 



CT+ I 



0<Z0 



<T 



t t 



(» + !)!! 



, (2cr— 1) ! ! 7T. 



Collecting the results, we may put generally 



f + ' 7 (» + <r-l):!(>-2)!! .„ , , . 

Q; dp = c<r 7---TT7;-,- , v , , it »i + <r be even, 

where c is equal to 2 or tt according as cr be even or odd, When (n + o) is odd, the 
integral vanishes. 

The results of this paragraph allow us to represent a quantity which is constant 
over a sphere, in terms of a series of tesseral functions which are all of the same type, 
the type being arbitrary. It follows that we may express cos crX and sin <j\ in terms 
of a series of tesseral harmonics of type cr. 

v on. co. "A. & \j 
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To find the values of the coefficients of such a series we put 

cos o-<f> = A£T* + a: + ,t: + , 4- • • • A£T; + 



a & a 



and integrating over unit sphere, after multiplying as usual by T£, we find 



2irr + l 



J -1 



Q^ cos 2 acf) d<f) dfi = A 



•2tt r + i 



(T;) 3 d<f> dfji = ttVt • 7 fi A - 

.o j -l 2?i + 1 (n — a)\ 



Substituting on the left-hand side the integral found above we obtain 



A" 



2n + 1 (?i 



co- 



JwJ 



<xj I (?fc + o- - 1) ! ! - 2) ! ! 2?& + 1 (-??, — <r — 1) ! ! (w — 2) ! t 

(n -j- a) ! (?t - a) ! ! (■/& -f 1) ! ! ~" 2 ' (?t + ff) ! ! (w + 1) ! ! ' 



where c is 2. or tt according: as <j and n are both even or both odd. 



-£> 



We obtain, for instance, in this way for cr = 2. 



o 



.J. 5 «j , -J . -fc 



9 



*j 



13 



■x - , « „ . <*■> t 3 . 4 . 5 . 6 •** "T" 5 . 6 . 7 . 8 ~ v 6 "*" • " ' 



+ 1 



J ~i 



Q* sin^ <//x, 



ftQn sin A l9 djji. 



If we put p = <T + X -f 2 in equation (G), and after changing from n + 1 to n 
integrate both sides, we obtain : 



+i 



_i 



sin* 6 Q* d[i 



n -j- a — - 1) (?i — X — 2 

(?i — Gr) (?i -f- X -f 1 



' + 1 



sin* Q;U ^ 



j „.j 



_ (n + or -- J.)_(>_+_5 - 3) (w - X - 2) (?i - X - 4) f + 1 

"" (w - <r) (a, - cr - 2) (n, + \ + 1) (n + X - 1) j _ 2 W4 ^ 

If n — cr be odd ? a continuation of this process reduces the degree of the tesseral 
function on the right-hand side until it becomes smaller than o\ and this will happen 
before any of the factors in the denominator become zero, so that in that case the 
integral on the left-hand side is zero. 

If n — <r be even, the integral on the right-hand side ultimately becomes 



f sin* 6 Ql dp - (2<r — 1) ! ! rsin< A+ ' +1 > d0 = (2<r - 1) ! ! v~^4!-,- c, 
where c is 2 or ?r according as X + a is even or odd. 



The integral to be determined now becomes 



r + i 



Th 



Ql sin A 9 dfi 
a be even) 



c - 



n + (T-l)!!((x T X 



(7i — o-) ! ! (?i + X + 1) ! ! 



(n 



2) (?i — X — 4) . . . (cr — X)] (if 



if n _ cr be odd. 
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A little care is necessary in the interpretation of the square bracket on the right- 
hand side. All the factors of the product contained in it may be positive or 
negative, but when X is intermediate between n and cr, some may be positive and some 
negative. In the latter case, one of the factors will be zero when or — X is even, and 
in that case the integral on the right-hand side is zero. The above expression does 
not include the special case n = cr. By extending the notation of double factorials 
to negative numbers as defined in § 1, we may also write, including the cases when 
n — X, or both n — X and <x — X, are negative, or when n = a, 

f+i . Kp 7 (*i+0--l)!!(er + \)!!(rc-X-2)!t 

L x Q * sm Ba * ^ c (» -^T^^^ lf n ^ a be even > 

= if n — • a be odd. 

r+i 
The integral fjiQlsm K cIjjl reduces to the one just determined with the help of 

equation (A). We thus find : 

+ 1 a, . A /)/ (?& + <r)!!(<r + X)!!(u~X-3)!! .. . ._ 

= if n — cr be even. 

The factor c takes, as before, the value 2 or tt according as a -f X is even or odd. 
For the special case a = 0, the tesseral harmonics reduce to zonal harmonics, and 
the last equation becomes 



f 



^ an tf ^ = c __ I ^__^^__ 1 ^ lf „ be odd. 



If X be even and n > X + 2, the fraction, on the right-hand side is zero, because in 
that case ( — X -— 2) ! ! = cc , and the numerator remains finite, but if n < X + 2, the 

value of -7 — — — 7--~ remains finite whether X be even or odd, and, in that case, the 



(-X 

right-hand side may also be written (avoiding negative arguments of the factorials) : 

, s^zl n ! ! X I ! X 1 1 



(n -1)11 (n + X + 2) ! ! (X - a, + 1) ! ! 



If n be odd and n s^ X + 2, we may transform the negative factorial and write the 
value of the integral 

— / — lY^ 1 tt itXl!X!i(tt--X--3) ii 

Similarly we obtain 

2 B 2 
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. „ sm A du = c -777 — — ^77T7~~ — "~77t;j it ^ be even, other 

«! r ?i !!( — X - 2) !!(« + X + 1) !! 



wise the integral is zero. 



This includes, as particular cases, 

+i 
P„ sin A 6 dfx ^ if X be even and ?* > X + 1 



r + l 
'-I 



( — 1) 2 7T ' -^ L -777--TTVT^r— "^ lf x be odd and ?* 5: \ + 1 

x ' n I ! (n -f X -f .1.) ! ! 



/ 1 \ -r \ 71/ — 1) ' ! X ! ! X ! ! , p x 

( — 1) - tf T777'~7 77777: ""77. lf : n < A + 1. 

v ; ?i ! ! (X + n + 1)!! (X — rc) ! ! 

Dr. W. I), NiveN" ( c Phil. Trans,/ vol 170 (1879 L), p. 379) lias already obtained 

r-hi 

an expression for the integral P n sin* 6 dp,. His results are identical with the 

above, when allowing for the difference in the notation and after correcting an obvious 
misprint in the equation marked (16) on p. 888. 



r+i ^pp 

If we write down the differential coefficients of P w . by means of the equation 

dV„jdp ~ (2n — I) P n ^.i + (2n — 5) P,,_. 3 + - • »? 

and repeat the same process p times, we obtain a series beginning with 

d^jdp? = (2ji - I) (2u - 3) . . . (2a ~ 2p + 1) P rt _ p + . . , 

There being no term containing a aonal harmonic of higher degree than n — p, we 
conclude that the above integral vanishes when e > n — p 5 and that when e = n — p 

' +1 o r?T->, 7 2(2n~l)l\ 

,, r? df^' U ^^(2n-2p + 1)1 ! 

If € < 7i — p we may transform the integral as follows : 



U € <v ^ ^ 2"^! L, " "V - * ^ ^ 

After € partial integrations, in which the first term always vanishes at both limits, 
the integral becomes 
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The integral on the right-hand side has been found in the last paragraph. Writing 
p -f e = a and e — p = X, we note that (n — X — 2) = (n — e + p — 2) is necessarily 
positive as e < w — p, and that cr — X is positive and even. If, further, w — cr be 
even n — \ must be even also, because the difference between these quantities is 
a — X = 2p, We may now write the value of the integral 



+ 1 jp p 

p. v^ ^ 



" rfuP" 



9 



7i 4- e -f p — 1) ! ! (vt — e 4- /o — 2) ! ! 1 



//, 



e - /o) ! ! (ft 4- e - p 4- 1) ! ! (2/o - 2) ! ! 



if 



e 



71 



0, in all other eases. 



If p = 1 ? the above reduces to 



+i 



rfP 



P e --— (//x = 2. 



p and n + p + e is even, 



§6. f +1 Qr-Qrc/fi, f +1 Q;Q;d/i, [ +1 QjQ^o^. 

•'-i j -i --I 

Before discussing the general integral of the product of two tesseral functions, it is 
convenient to obtain the solution in a few special cases. When the type of the two 
tesserals differs by 2, we may transform the integrals as follows : 



+i 



(i 



<r+2 ,73- + 2 p <r r iv p 



dftf 
+ 1 



tf/ 



+ 1 



<T 



-ft 



<r + T M rn\- 



.7 (T + 2 



dp, 



9-, 



i 4- 1 J 



(i - ^) , 

•1 Clfl 



dfjb 

2V+-1 tt -1 n I "• ± . ( - + i 



rf 



o* + i 



\ <i/x ( 






c//x. 



By the application of Rodriguez's theorem the right-hand side reduces to 



(-l)' +1 '(i + o-)!f +1 d T+2 P M 



rf' 



'CT— I"!) 



(i +<r +2)(i+(r+l)'i^§±i 



(i — cr) (*'— cr— l) 



cl 



a— lp 



P 

-M-l 



6?/^ 



or+1 



cfyx. 



Integrating each term partially or + I times we arrive at the equation 



+i 



y w " ^f ^M 



1 (i 4- a) 



+UIV 



n [(i + a + 2) (t -f cr + 1) P m - (i - cr) (i - cr - 1) P^J ^ 



2^ 4- 1 (i — cr) ! * -i ^/x 

If i + ft be odd, or if i > n y the integrals on the right-hand side vanish* 
and i -\- n even, we must substitute, 



If i < n 



r+irfp M C +l dP n 

~ Phi = T^ Pi-i = 2, as found above, 
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which gives 



+ 1 



Qr 2 Qr^ 



4 -; , J r~ { (T -f" 1 ) # 

% — a) ! J 



When n = i, that part of the integral which depends on P i+1 vanishes, so that in 



that case 



+i 



Q^QZ dfi = 



('?l + cr) ! 



2n + 1- (> 



<j 



2H 



ore particularly when <x = we have the integrals 



+i 



QlPidfi = 4, if n -f- i is even and i > n 



in all other cases, 



r+l 



-l 



*%«. -t «W'^, 



2- 



n . 7i 



2% + I. 



The result that 



+i 



Qn +2 Qidfi has zero value whenever i > n can be extended to 



the more general integral 
only consider the series, 



+i 



J~i 



QnQUfy- provided then cr > p. To show this we need 



Ql ~— A Q« 2 + A 2 Q^ + 2 4" 



& » * 



where the coefficients may be determined from the integrals found above. Multiplying 
both sides with Qf" 4 and integrating, all the terms vanish when i > n and hence 

Qn Qf~ 4 = in that case. From or — 4 we may proceed to cr — 6 and so on. 

To obtain the integral Q^Qjcfyt, we llse formula (F), multiplying both sides by 

Q p n and integrating. We may without loss of generality take p to be smaller than a. 
Utilising the result which has just been obtained, one integral on each side drops 
out and the equation becomes for even values of or -f- p, 



■+i 



Q«Q'° 



•+i 



( n — a- -±. 1} ( n _ a -j_ 2) j Q£ 2 Q^//x 



(n — cr + 1) (n — cr + 2) (n 



(j 



n 



cr + 4) | U Qr 4 Q£4* 



(-1) 



( — X ) 



<r- 



(u — o-) ! . 



£=£ (?i + p) ! 



Qftfodp 



9, 



{n — cr) ! " 271 + 1 



if p < <x and p + or even 



if p -f- a* is odd. 
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Reverting now to equation F and multiplying both sides with Q£_ 2 > we obtain 



+i 



Wi / nS c in—2^ f^ 



( n _ a- + i) ( n _ „ + 2 ) f + Qr 2 Q^ 2 ^/x 



+ 1 



+ (n + cr - 2) (n -her — 3) j Q^Q^p, + f Q*UQUdu 



• + 1 



(w ~ o- + 1) (n - o- + 2) Qr 2 Q p B -s^ 



cr — p + 2 »,) 

+ (-1) s , 



76 — o 



('/i + cr —2) (n -j- cr — 3) - 



7i — 2 -J- p)! 

(n — a) ! 



<> 



n ~2 + p) !" 
<> - cr" 2) ! 



r+i cr-p+2 / w i „ 9\ 

= _ ( n _ o- + 1) ( n - o- + 2) j_ Qr 2 Q^- 2 ^ - (- I) * 4 (cr - 1) l • J^- , ,~ j 

The integral on the right-hand side may again be transformed in the same 
manner, changing from a to or — 2. This leads to 



f +1 Q;Qf,-2^ = (n-cr+l) (n-a + 2) (/i-cr + 3) ( 



?2 



cr+4) f +J Qr 4 QU^ 
--1 



or — p + 2 // 1? i n »A | 

+ (-1) * 4 ( ''^^ ) ; ) '{(ct-1) + (o--3)}. 



If the same process be continued until the integral on the right-hand side 

r+l 

becomes Qu +2 Q^~2^/ X 5 the remaining terms on that side will consist of a series in 

arithmetical progression. 
Adding this, we find 



■+i 



»i 



"Q;Qi-^ = (- i) 



<T~p 



+ ° 



"('«• 



9 



P)' 



(?6 — a) I 



+ 1 



Qf 2 Q>w^ 






4~ 2) (cr — p — 2) 



o—p + 2 
(-l)"2~ 



(n - 2 - />) ! fa - 2 4- /?) ! 

(?i — <x) ! ' (76 — 2 — p) I 



(4p + 4) 



6^ — 2 -f- s ) I 

_j ^__^ ' (or + p 2) (cr 

(n — a) ! x ' r / v 



9 



) 



(— 1) 2 — •-■, ------ t ••:— ' (cr 2 — p 2 ), if cr + p be even and p < cr. 

\ / {^)i _ ^r) , ^ ' 



If o- + p is odd or if p > cr the integral is zero. 



192 PKOFESSOJft A. SCHUSTEK ON SOME DEFINITE INTEGRALS, 



§7. fQsP.-tfyG, f Q;Q?rf/x. 

The integral of the product of a tessera! function and of a power of sin 6 having 
been obtained in § 4, the above integrals are found by expanding the tesseral 
function or the zonal harmonic in a series proceeding by powers of sin 9 

From the expression of a tesseral function of degree n and type or, as it is generally 
given 9 viz, : 

sin tf [A /x ;t ^ + A^'*"' 1 + , . . ] s 

it is seen by writing x = sin # ? /x = \/l — as 3 , and expanding the binomials, that the 
term of lowest power will be sin " > and that if n — cr be even, the term of highest 
power is sin w $. 

The coefficients of the series, which are given by Heine, are most easily obtained 
by going back to the differential equation : 

d ,. „ d ' 



Substituting x = v 1 — /x 3 , and changing variables, this becomes 



n . ra + 1 . Q; + - Vl1 + ( 1. — ^) - 7 V- = o Q'. 



If the series 



a^ " + a <r+2 a; ' +3 + • • • «,,# 7 + a, /+i ,a:? + ~ + . . . 



satisfies this differential equation, the coefficients a qvti and a f/ must be connected by 
the relation 

a q (n — q) (n + # + 1) + a, /+2 (y — cr + 2) (q + cr + 2) = 0, 

as is seen bv substitution. 

The first coefficient is determined by the fact that it is equal to the value of 



Hi "4™ (T ) ' I 

d'VJdu? when a = 1. This quantity is known to be equal to ,-'■ *ir-\" . The 

other coefficients mav now be determined in terms of this, and we find in this way 
lor vvj^ une soi ies 



O - cr)I ( 2o-)f! 






1 cr + 1 V 2 / 1.2 (7 + 1 . a + ^ \ ^ 



The series breaks off with the term x n ~*> when n — or is an even number, but also 
holds if n — cr is odd. 

The factor of x k in the series reduces to 
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1 ' } r n -a- !)!!(> -X)!! (X - a)\\ (X -f a-)!! 



When (u — <r) is an odd number it will bo more convenient to use a different 

series. Writing Q = /xN, and changing the variable to x = \/l — * /x 2 as before, the 
differential equation becomes 

n — 1 . n + 2 N v ; -7- + (1 — re 2 ) — - = -- - , 

x dx v y aar or 

and from this we find for Ql the series 

( ?l ^ (J )!(2o-)!!L I. crH-"r V-V ~"~" ll 7+T/7+2 \£ 

The factor of fix k in the series is found to be 

\*(n + <r- !)!!(?& + X)!! 1 



(-1) 



71 - a-)!! (r/.. - X - 1)!! (X - a)!! (X + a)\ 



r + 1 

In considering the integral QjlQ? fifyt, we may, without loss of generality, take p 

* — 1 

to be smaller than o\ If of the two quantities i — p and n — <x, one is odd and the 
other even, the integral vanishes, 

If (n — cr) and (i — p) are both even numbers, we may express Q* in terms of a 
series of powers of sin and obtain in this way : 

- S-W- 1 , A :7(^ + ^)n(^4- X~-l)!I (74-p ~l)!! (7-X^ 2) ll( p + X) !I 1 



A = cr 



The symbol 2 (J) is intended to express that X takes successively the values 

A -scr 

or, cr -f 2, . . . 7?, leaving out the alternate numbers. The constant c is equal to 2 or tt, 
according as p — cr is even or odd. 

If 7i — cr and i — p are both odd, we find in the same way : 

+ 1 -. A~n~l *-"<r()l 4- cr — IVtfe 4- XV t I f + i 

— 1 A = - V^ — <Tj : ; ^<?£ ~- A. — J. J : ; {A.- — (jj : ; (A- -f- CTJ j : J _ 1 

*^y i^^ 1 - 

The value of c is the same as before. 

Our result shows that when p — cr is even, the integral vanishes when i > n. 
For when p — <r is even, p — X — 2 will also be even, for all values of X, and 
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It will be negative because or is the smallest value of X and a is larger than p. 
l/(p — ■ X — 2) ! ! must therefore be zero, and all the terms of the sum vanish unless 
there is a negative even factorial in the numerator. The only factorial which can 
become negative is (i — - X — 2) ! !, but n being the highest value of X, it will not be 
negative when i > n. For the special case, % = n> all the terms of the sum vanish 
except the last, for which X = n> We have in that case 



(i- 


— \ «— 2) 1 : 


(i ~~n- 


-2)11 _ (-2)!! 


(p- 


— X -— 2) ! i 


(p — n - 


- 2) ! ! — (p - rc - 2) ! ! 



n—p 



| rasa J. S " I 11 tsaasas * P ) i 



^-TiL 2 (% 4- p) ' 

and the integral reduces to {— 1) 2 -- ••-•-; -\-f, as already found in 8 6, 



The same expression is found when n — cr is odd. 

Another and generally more convenient expression f 

expanding Q? instead of Qj in terms of powers of sin 
n — cr and % — p are both even, 



or QnQ?^ is found b 
•'-1 

$ s We thus obtain, when 



1 4- 



X=i 



w>,7 _.^/ n -~ p -(» + ^-l)"(*L+/i>)n(i + X-l)!!(»-X-2)!! 1 («r + X)!! 

l ^i'"/*- 6 - I l ) - ( n _ ff )!!(i-p_I)!!(t_\J!!(„ + \ + l)!!(X-p)!!(\ + p)!!( < r-X-2)!!' 



and when n ~ a- and i — p are both odd, 



A=i-1 

c 2 (2> ( 



ritmg 



1) 







__ (<7 + X)!! 



^"e(7i + ^)!!(i + /o-l)!!(^-X,-3)!!(i + X)!! 1 



(cr 4 p) (o- 4- p — 2) , . . (a* — p 4 



2) (cr — p) 



(2p)! 



O 



(cr 4 p 4 2) (cr 4- p) . . . (cr ~~ p) (cr — p — 2) 

(cr + p 4 4) (cr 4 p 4- 2) . . . (cr — p — 2) (cr — p — 4) 

411(2/3 4- 4)!! 



we may put the integral into the following form 



where 







r+l 

• -1 



iQ^dii, if (n — cr) even, (i — p) even, and excluding n = i 



- c 



i 4 p) 1 (% 4 cr — 1) ! 1 (w 



<>\ » t 



(i - p) ! (n - <x) ! ! (n 4 p 4 1) 1 



\ f t ^j 



Ac 



(i + p + 1) (i — p) 



- f % __ p — 2) (?t 4 p 4 3) 



A. 



i 4 p 4 1) (i 4 p 4 3) (i — p) (ft — p — 2) 
7l , _ p — 2) (w — p — 4) (> -I- p 4 3) (w 4 p 4 5) 



wnere 



(3 



= 2, if cr — p is even, and c = 75% if cr — p is odd 
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the last term of the series is 



[(i + p + l)(i + p + 3)... (2i - 1)] [(i - P ) (i - p - 2) . . . 2] 



*-/> 



[(ft — p — 2) (?i — p — 4) . . . (» — i)] [(ft + p 4- 3) (ft 4 p 4- 5) . . . (w 4- i + 1)] ? 



but the series breaks off before the end, if cr — p is even and i > cr owing to the 

factor A vanishing. The number of terms in that case is cr — p — 1, 
Similarly when (n — cr) and (i — p) are both odd 



•+i 



^%nSti ® / l L 



C 



(i 4- p) ! (» + cr) 1 1 (ft — p — 3) ! ! ^ 

(i — p) • (ft — cr — 1) ! I (ft 4- p 4- 2) ! ! """ 



wnere 










A ft 4-p-r 2) (z — p— 1) . (i4-p4- 2) (i + p + 4)ft~~ p — 1) (i — p~~~ 3) 

2 (ft — - p — 3) (ft 4- p 4- 4) '* (ft — p — 3) (w. — p — 5) (ft 4- p 4- 4) (ft 4- p 4- 6) ' * 



where c and the coefficients A have the same value as before. The last term of the 
series is 



__Kj ±P ± ,2) (i + p + 4)...(2j - 1)] f( ; i + /) + l)(i«.p- 3) . . . 2|__ 

l ^ J [(ft - p - 3) (ft - p - 5) f. . (^ - "*)]'[(» + P T4)7»TpT6y.T/(»"+~i +T)J 






"fo obtain JQ£P t - d^ we need only put p — in the previous investigation. This 



gives : 



A «— « C/ » 



cr + 2 . €7 . cr ~~ 2 



A 



o~ 4- 4 . cr 4- 2 . 



cr . cr 



. (J 



Mm 



4 



M » ^C Hi iLJ * -£ 



f 



+1 f>p 7 _ 9 (^ + cr-l)!l(ft-2)!l T t.j + 1 |A "i.i-2.i4-'l.i4-3 

-l (n-«r)!!(»+l;!! L 8 »-2.« + 3^ *n-2.n-4.ra + 3.» + 5 



■f> » * » i «( 



if * and cr is even, and 



■i * ' p (ft-cr-l)!!(/i4-2)!I 



A 



o 



A i — 1.^4-2 . . ■£ — l.i — 3 i~{~ 2.^4-4 



-2 



ft — 3 ft 4- 4 



7& — 3.72, — 5 ft 4- 4, ft 4- 8 # * * 



if or is even and i odd. Thus, if i be even 



'+ ,i. 



O 2 P* diL 



4 



a. 



Q/P f - cZ-u, = 8 [(n + 3 . n — 2) — 3 (i . i + 1)] 



■+i 



-i 



gT> 



P 2 - % = 12 [(n + 5. ra + 3. n — 2. n 



■m -*- z 



4) — 8 (i . i + 1 , n + 5 . n — 

+ 10 ({ . i — 2 . i + 1 . i + 3)] ; 



and if i be odd : 
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r + l 



'+1 

■1 

'+1 



^L i 



Q/P, (J [i = 8 [(ti + 4 . n — 3) — 3 (i + 2 . i — 1)' 



Q/Pi tfy, = 12 [(n + 6 . n + 4 . n 



+ 10 (i + 2 . i + 4 . i — 1 . i - 






The equations do not hold for i = n as has already been explained. 



8. I Q* sin p0 d!/x, 



Q£ cos p0 (Ifi, 



These important integrals are obtained in two different forms according as Q£ or 
the trigonometrical functions are expressed in terms of the powers of sin 0. In the 
former case the problem reduces itself to the evaluation of integrals of the form 



7T 



7T 



sin x 6 sin pd dp and sin x cos # sin jp# c?0. 

*' ■» 



It may easily be proved that 



+ 1 



p-l 



sin x i9sin j3#c//x = (— 1) 2 c 



(X 4* 1) ' 



(X -f 1 + _p) ! ! (X. + 1 



- — . when p Is odd, 
— p) ! ! 



when, p is even, 



ii 



sin* cos p0 rZ/x = ( — 1 



(X 4- 1) ! 

G (x"+ TTp) TT(xT" i 



~ :-,-; when JO Is e 



when p is odd, 



sin*" 1 £ cos $ sin p0 e//x = ( - 1 f*~cp ■ ^ ] + v) , , ~ + y _ - p) ; , when p is even, 



when p is odd, 



+i 



sin x 1 cos cos p# ti/x 



-p-i 



1) * cp 



x 



(X + 1 +p)\\(\ 4- 



when ?9 is odd, 

p) ! ! l 



when p Is even. 



In these equations c is equal to 2 or 7r according as (p + X) is even or odd. 
From the results of § 7 we may now write, if n — a* be even and p odd, 



f 
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+1 
QlsinpO d[i 

vr>/ ,y^(n + <r)\\(n + \-l)l\ 1 f + 1 . A/} . . 7 

= 2 W (— I) - ;— -77, TT7, /'; r,T7T TT, \ sm " Sill pV dil 

a^ct X ; {n~(T - .1) ! ! (71 - X) ! ! (X - cr) ! ! (X + cr) ! ! J^ z 

^ , v x+ ?:^=1 (n-l- <r)i!( tt + \-l)i ! (\ + 1)1 _ _1_ 

^ ^ } " (n-<r-l)!!(tt-\)r!(\-<r)!l(X + ^^ 

When p and n — cr are both even the Integral vanishes. 
If (n — cr) be odd and p even, 

+i 
Ql sin p0 clfjL 
■i 

= 2 W (— 1) 2 r r7T7~~~T~~Trrv /"T T77~T~~ — rrA cos v but sin pvdu 

A=(r x ; (n — a) Win — X — 1) ! ! (X — <x) ! ! (X -f cr) ! ! J _ x l r 

-™vV/ ^r^(n + a-l)\l(n + \)l\ (\ + 1) ! 1 

* A ~, l " ; (n^a)!!(^^X^l)!!(X^ -)!!(X + o-)!!(X-f2^)^(X+2---p)!! , 

It will be noticed that the integral is zero whenever (n + <r + p) is even. Hence if 
sin p# is expressed in terms of tesseral functions, p being even, only odd values of n 
occur when cr is even, and only even values when a is an odd number. The reverse 
will be the case when p is odd. 

We find similarly, if (n — cr) is odd and p odd, 

+i 

QjJ cos }) cIjjl 
-1 

= z ( - ) (— 1) - . -77 t-".^- ,a ^rrTTT^Trt cosffsm A 0oosp0duL 



° P ?- 1 ^ X ' " 0^"<x) ! ! 7^"-X JlTi"! fX-o-) ! ! (X + 



^ v ' fa - a) ! ! fa -X-l) ! ! (X-cr) ! ! (X + <r) ! ! (X + 2-f p) ! ! (X+ 2~p) ! ' ' 

and when ?i — cr and p are both even, 

Ql cos j? @ dfj* 

•i 

- . VS /_ ^^" (?l + ^ * * (^ + X«™ 1) ! ! (X -f 1) ! 1 , 

^^^ ^ ^ * ( ? i- o-- 1) ! !(w-X) ! ! (X- o-) ! ! (X + o-) ! ! (X+ 1 +.p) ! ! (X + l-p) M J 

when 7i — <r + P is 0( ld the integrals vanish for all values of p. 

We may now prove the theorem which has been mentioned in the first paragraph. 

- r+l 

If cr and|) are both even, Ql sin p6 dp is zero, unless n is odd. In that case the 
quantity (X + 2 — ■ p) ! ! which occurs in the denominator will be even ; as X takes up 
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successively the values cr ? cr + 2, &c. The highest value which X + 2 — p can take 
is n ~{- 1 — _p, and this is negative when p > w -f- 1. 



'+i 



follows that 



£ sin jo$ d/x = when p > n + 1 and cr is even. 



The restriction that p is even is not necessary, because when p is odd n will be even 
and the factorial in the denominator, according to the above equations, is (X-f-1"— p)\ ! 
This again will be infinitely large whenever it is negative, or taking the highest 
value of X, which is now n, whenever p> n + 1. 

We prove exactly in the same way that 

Ql cos p0 djjL = when p > n + 1 and cr is odd. 



The value of the integrals in the general case may conveniently be expressed by 
two series Si and So defined by 



(ft - p + 1) ! ! (ft + J? -f 1) ! ! (71 - cr) ! ^ A 






(71 -f (r) (ft — a-) (ft 4- p + 1) (ft — jo + 1) 



2% - 1) ! ! (n + 1) ! - ,J - — * 2 . 71 + 1 . ft . 2ft •- 1 

(n 4- c) (n -f o- ■— 2) (n — cr) (ft — cr — 2) (w 4- p 4- 1) (ft 4- p — 1) (ft — p 4- 1) (ft — £> — 1) 



2.4 



ft 4- 1 . ft- ft — 1 - ft — 2 



_/?'/£ — _£ ^ *^7t^ "**^ O 



8 © S 



(ft — p 4- 1) ! ! (ft 4- _p 4- 1) ! 1 (ft -^<r) ! ^ _ . _ (^ + ^— 1) (^ — °" t 1) ( w _+.^ + *) (^ — P + 1) 
(2ft — 1) ! ! n ! j? "^ 2 . ft . ft — 1 . 2?& — 1 

(ft-f cr — 1) (ft4-cr — 3) (n — a—l) (ft — cr — 3) (?7, 4-/>4-l) (ft+jfl — 1) (ft— ^? 4-1) (ft— p — 1) 






9 



% , n 



n — 



9 



ft 



O 



^'ft ***""'* X *dft ~™*~ O 



* « 9 



T1 
her 



+ 1 n '~°" + ^"::i 

!£ sin j;$ dfi = (■— 1) " c 2! if p be odd. and n + cr is even, 
■1 



(- 1) 



n — tr+'p-Y 1 






; So if v be even and w + cr is odd. 



whenever ?i + p + cr is even, 



+1 



Q£ cos jp# fi/x 



72- — (T + 'j) 



1 ^ 



c S] if p be even and n + cr is even. 



'/l — (T + 7>+2 



( ~*^ 1 ) 



c S 3 if p be odd and n + cr is odd. 



,u rn a. 



wnenever n -4- p ~j- cr is oaa. 

The constant c takes the value 2 or tt, according as p + r ^ is even or odd. The 
series. Si and So. break off in all cases ultimately when one of the factors of the 
products (n = -cr) .. (n — cr — 2) , » . or (n — cr — l) (?i — a* — 3) . . . become zero. 
When n — p is odd, the series may break off before it has reached its full number of 



WITH APPLICATION TO SPHERICAL HAEMONIC ANALYSIS. 



199 



terms, owing to one of the factors In the product (n — p + 1) (n — p -— 1) , . 
becoming zero. In the special ease n == a the value of l£ l reduces to 

2n - 1) ! ! (n + 1) ! 



7i — p 4- 1) ! ! (n + p ~f 1) ! ! 



The above equations are inconvenient for numerical calculations, unless n ■— <x be 
small, or n -— p be even and small. We obtain quite different and generally more 
convenient expressions, if we begin by expressing sin pO or cos _p# in terms of a series 
proceeding by powers of sin a 

Let 






X . jU 



Jl . ^w . tJ 



and generally 



B 



A 



O - \ 4- 1) (p - X -f 3) , . . (p + X -- 3) (p + X — 1) 



X! 



(p + x - 1) ! ! 1 

0* - x - i)Ti xl 



Also put 







" ^ 1 ; (\ = p . ; Co = p . | ; C 3 = p 



jo-l .p+ 1. p _ p-2.p.p + 2 



1 9 '■>> 



4! 



C, 



79. 



P-X + 2) (p — X + 4) . . . Q? + X — 4) (^ + X — 2) 



X! 



~ xi' '^(7^)Tr 

Then the well-known expressions for the trigonometrical functions of the multiples 
of an angle may be written : 
If p be even : 

sm p - — B, sin $ — B, sin 3 6 4~ B, sin 5 6 - 

cos l ° ,} 

cos p$ = C — Co sin 2 + Qj, sin 4 ' $ 
and if p be odd : 

sin p9 = Oj sin — C<> sin 3 + C 5 sin 5 - 



o » » 



db Vi sm 7 '""" 1 0, 



, . . i C„ sinJ° # ; 



± CL sin' 0, 



cos pO 

COS 







B 3 sin 2 6 + B 4 sin* 



. . . ± B^jl sin*"" 1 6. 



We derive from this, p being even : 



-i 
+ i 






A-l 



QU sin pd dfi = 2 (2) (- 1) * B A j Q£ sin x cos dp, 



\=l 



■ + 1 
•1 



Q; cos pd dfji = $&■ (- I) 2 C A f +1 Q; sin A d^ 
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and p being odd : 



'+1 



Ql sin p9 dp 



A = l 



A-l 



l) a c 



+ J 
A ! Q£ sin A 9 dji, 



4-1 



Ql cos p6 dfi 



\zzp~l 

. V(2) / 

A=G 



A 



+ 1 



1 ) B A Qf h sln A $ cos $ dfi. 



Hence from the results of § 4, if jp be even and (n — or) odd 



_i * r r a=i ; — cr— 1)!! (w + A, + 2)!!(o- — X — 2)!! 



The constant c is 2 or b\ according 1 as <x is odd or even 



fc> 



If p be odd and n — er even, and with the same meaning of c : 
Similarly if p and n — cr are both even, 



, +i 



Q^ COS p$ dfi 



c %w ( 

A=0 



A 



1 (ti + <y — 1) ! ! (?i — X — 2) ! ! (or -f X) ! ? 



and if p and n — ■ <r be both odd ? 

Q* cos pOd/A = c S (2) (— 1) *" 

-1 ' ' A-0 



(71 + cr) ! 1 (n — X — 3) ! I (<x 4- X) ! ! 



A (^ 



(X 



1) S ! (n + X + 2) ! I (<r - X - 2) ! ! 



In the two last equations c is equal to 2 when cr is even and equal to it when a is 
odd. In the interpretation of the summation a little care is required when cr— X — 2 
and ?i — X — 2 or n — X — 8 are negative, as the alternate factorials are infinite for 
even values of the argument, 

To put the equations into a form useful for numerical calculations, it is convenient 
to designate by separate letters the following six series :— 



* cr "— 1 . cr -f- 1 "•- C a 



+ 



er 



n 

O . (7 



1 . cr -f 1 . cr + I 






n 



»> 



*T* *x 



or -— o . cr 






cr 



l.cr + 1.0* + 3.er + 5 



71 



- 3 . n — 5 , n + 4 . ?i + 



<j' 



I\ * cr — 1 « a 5 4* 1 *"" 1^ 



J_ 



o 



'B 5 ^ 



cr 
n 

cr *— 5 . cr 



1 . cr -f 1 . cr + 3 
4 . % -f 5 






CT 



l.cr + l.tr + S.^-fo 



n 



4, 



t~7 



I- . % — 6 . n -f 5 . % ~f / 



• » * 






O w 



a 



cr — 2.cr.cr+2 . , N cr — 4 s 0" — 2 , cr » cr -f 2 , cr 4- 4 

; jL i> _ - 

''" n — 2 . n + 3 ' 4 ^ — 2 . ti — 4 . 71 + 3 . 71 + 5 
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y «r = T> _ T> 



0' 



cr — • 2 . or . cr + 2 jy a — 4 . cr — 2 . cr . cr 4~ 2 , cr 4~ 4 

7i — 3 . w + 4 4 % — 3 . % — 5 . w -f 4 . n. + 6 



« » # 






7i 4~ 2 — a- . n 4- 2 4- cr 



#+2 



J , 2?Z -f* & 



+ 0^ + . 



?i + 2 — a . n 4- 4 — cr . n 4 2 4 o" . ti 4- 4 -f o* 

2,4 2?^ 4 3 . 2?i 4- 5 "™ 



N> rz: "R _ "R 



n + 1 — cr , ?t 4- 1 4- cr 

.4 • <y?t 4~ O 

71 4* 1 — * 



o 



B 



cr . % 4- o — cr 



rt+3 



, n 4- l + ^.w-fS-ftr 



2.4 2?i 4 3, 2?& 4 5 



« * « 



The series are all continued until they break off, which may happen either 
the factors B, 0, or one of the other factors takes zero value. 

With the help of the functions just defined, the integrals now take the form 



se 



+i 



Ql sin 2^0 d[i ™ rr 






(n 4- or ~ 1) ! ! (yl - 
(^ - ex) ! ! (71 4- 2) ! ! 



M* if tr be even, p odd, n even, 



(n 4- <r) : i O — 4) ; i Q . ,» , 

^ __ r-oCTT few " °* be even, w even, n odd, 

(n — cr — 1) ! ! (pi 4- o) ! ! ' £ 



o 



'/& 4- cr — 1) ! ! (n — 3) ! ! »rJ ( w . 4 cr) ; Tr n 

?6 — CTj i : {'ib 4 J J ; : (^?/, 4 1 j » ; 

if cr he odd, p odd. 71 odd. 



+1 



2 



r?i 



> 4- gQ ! ! (71 -4)!! , 

~cr l)!!(^ + 3)'!! ' l h ^ 



if n -f* cr + _p be even, 



' (2b + 1)!! w 

if cr be odd, p even, n even. 



Q^ cos p0 c h L 






C'l + o-_-l)_!|(«_-_2) ! ! T> , . :- (n + <r) ! T<v 

( n _ o-) : : ( re _ 1 ) ! ! L ' " "♦" ^ ■ / (2n + 1) ! ! " 






n 4. o-)!!^-, ,)!! 



(a — cr 



1) ! ! (u 4 2) ! 



T^r; V; + ( 



if cr be even, ^ even, n even. 

' ; (27i 4-1)!! ^' 
if a be even, p odd, n odd. 
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tt^-— — -— — — ^--r— - — UJ it cr be odd. p even, 7Lodd. 

(n — a) 1 1 (n 4 1) ! S x 

(% + (r) ! ! (?i - 3) ! ! tt • p i Ti it 

^^ _ -— - V^ it cr be odd, p odd, n even, 

(71, — cr — 1) ! ! (??, 4- 2) ! ! n ? r 5 

if n + o- + p be odd. 

2 D 
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By comparing these last equations with the results obtained when the same 
integrals were expressed in terms of the series 2q and S 2 , we may derive the following 
special values : 

M£ = if a be even, p odd, n even and p = n + 3, 

S^ = „ cr „ ,, jp even, n odd „ p ~ n + 3, 

U^ = „ a „ odd, p „ n „ „ j^g^ + 3,, 



V^ = cr 

y ,>. w ?? w ?? ?? 



£> odd, ?? even ,, jp ^ n + 3. 



The two series marked N* and K* had to be introduced because the series 
for M*, S^, U*, V* break off as soon as one of the factors becomes zero. But in the 
original summation which gave rise to, e.g., M^, viz. : 

A £|> (n -f <r - 1)1! (7i - X - 2)1! (a + \)!1 
^ ^ A (ti - o-)lT(?i + X 4- 1)!! (7- X -2)!! ' 

(cr — X — 2)!! begins to be infinite when X = cr, and hence the terms of the series will 
drop out until (n — X — 2) is negative, ?'.c\, until X~7?. For higher values of X, 

.--■ ^-~ will be finite again. There is, therefore, a second portion of the series 

(cr — \ — 2)!! 

not included under M^, and it is this second portion which appears as K* 



§ 9. Relations between M*, S*, U£, V* 
Certain relations exist between the series which serve to express the integrals 

+ 1 r + 1 

Q^ cos p# dfx and Q* sin £># c//x, and these relations are important for the 

calculation of numerical tables. 

Starting from the equation (F, § 2) 

Q* - Qj-2 = (ri + o- - 2) (n + cr - 3) Q*:l - (ri - cr + 2) (rc - cr + 1) Qr\ 

we may multiply by cos pd or sin pd and integrate on both sides. Expressing the 
result by means of the above series, we obtain the following set of equations : 

(n + cr - 1) (tt — 3) M; — (n — cr) (n + 2) M;_ 2 

= (n + cr - 2) (w + 2) M^i; - (n, — cr + l) (n — 3) Mr' 2 



(n + cr) (n - 4) S* - (n 

= (n + cr ■ 

(r, -fcr- l)(n- 2)U; 

. = (n + cr ■ 



cr — l)(n + 3)S^_ 2 
3) (n + 3) S;i2 - (n - cr + 2) (n 



4)S 



o— 2 

71 



(n-o-)(w + l)U;_ 2 

2) (n + 1) U;:i - (n - a + 1) (n - 2) Ur 2 



(n -j- cr) (n 



3)V^-(n-cr- l)(w + 2)V;_, 

= (n --f q- - 3) (n + 2) V^ - (n - a + 2) (n - 3) VjT 2 ) 



(K). 



WITH APPLICATION TO SPHERICAL HARMONIC ANALYSIS. 



203 



Other relations are obtained as follows : 



■+i 



IT 



TT 



Ql sin pd d[i = Qn s i n $ s i n P@ dd 

Jo 



d 



£> Jo 



cos p# ■-— Ql sin cZ$ 



d^ 



i r +1 d 

cos p# 7 Q* sin <ia, 



and making use of (H 2 )> § 2. 



1 r + l 

2p<x J _ j 



cos jq$ [(cr + 1) (n + cr) (n — cr + 1) Q'J * — (c — 1) Q» +1 ] d/x- 



If cr be even and n odd, the integral on the left-hand side depends on S£, while 
that on the right-hand side depends on U*"" 1 and U£ +1 . For or even and n even, the 
left-hand side may be expressed by M£, and the right-hand side by Vl" 1 and V* +1 . 

Treating the integral j Ql cos pOdfi in the same way, and collecting the results, we 

find : 



2po-Ml 

2pcrS: 



2pcrU; 



7i 



(cr + 1) (^ + cr) (^ — cr + 1) Vn l ~~~ (or — l) (n + cr + 1) (n — cr) V* +1 ] 

(w - 2) (w + 3) [(cr + i) ur 1 - (cr - i) ur 1 ] 

__ 1 

(n - 



o 



)(« + 3) K°" " !) ( ft + o- + 1) (« 



(or + l) (?Z + cr) (?l 



.j 



-|— 



i)sr 



2pcrV* = (cr — 1) 



cr+1 

"A 



(cr+ i)Mr x 



We may also connect together the different values which the same function assumes 
for different values of p. 
Starting from the identity 



■+i 



Ql [sin (p + 1)0-— sin (/? — !) 0] cfy, 



r+i 



9 



Q£ sin p6 cos d//, 



we find, transforming the right-hand side with the help of equation (A), § 2, the 



integral to be 



r + ] 



!;/, + 1 J 



■( n — cr + 1) Q* +1 + (n + cr) Q^J sinp$ rf/x 



Three other equations may be obtained corresponding to this ; by substituting the + 
for the — sign in the left-hand sign, or by changing the sine function into the cosine 
function. Each of these three equations gives two relations according as a is even or 
odd. The following eight relations are thus obtained : 

A~J XJ J-J 
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'~± (SJ {p + 1) + S* (p - 1)) = (^ - 2) M'l,, (p) + (w + 3) W-i (p), 

— — (M*(p + i) + Mu(jp - l)) = - r ii ^-- l H^ +1 (p) 



M 4~ 4 



+ -i^K^ g, ( }) 

""' + ] ' M r, (p + 1) - m^ (p - i)) = (ii + o- + i) (n - i) u:t\ (p) 

(,i _ ^ ( n + 2 ) Uji! (p) : 
2 " t ] ( S « (P + l ) "~ ^< (P ~~ *)) = ( w + °" + 2 ) ( H ~" 2 ) V '+S (?) 



(//, — a- — I) (n -j~ •->) V'i{ (p), 



2 - w .7 ■- (^ (p + 1) + u: (p - i)) = ( " + ff+ p^ --* + J) v; +1 (p) 



yt ~|- -> 



9 V '*-l UV: 



+ 9l _ 9 



2v ' + T (V; (p + 1) + Vj(p - 1)) = (n - 1) U:,, (p) + (« + 2) U^ (p), 



v>. 



^f x (Uj( P - i) - u:(p+ i)) = --^-J AM -i(p) - M-^M-Ky), 



§ 10. Numerical Calculations of the Four Series. 

For the calculation of the numerical values of the series M*(p), the special case 
a = was taken as starting-point. In this particular case the well-known integral 

which holds when n + 1> is odd and n ^jtf — 1, gives 

n ! ! (^ + 2} ! ! (n + p — 2) ! ! (71 — « — 2) ! ! , . 1 -1 , 

mo — _ • : ^!LZLiZ^_- ^ V— i-i ; - v *— — f— when 72 is even and « odd. 

The values of M° were calculated by this equation for all values of p up to and 
including jp = 11, and for all values of n up to and including 12. The first of the 
equations K, § 9, was then applied, and substituting a = 2, n = 2, and Mj; ~ 0, the 
value of M was calculated. Similarly putting n = 4, G, &c, the values of M; for even 
values of n were all tabulated. By successive steps I similarly found M^ ? M^ &c. ? up 
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to MJi The last value so obtained, and some of the intermediate ones, were calculated 
independently from the series defining M, numerical errors being easily detected in 
this way. The values of S;' ? IJj, V* were dealt with in a similar manner. As a 
final check, the values of S* and V* were obtained from U* and M* by means of the 
second and fourth of the equations (L). 



§ 11. Application of the Previous Results to the Expression of Functions by Means 

of a Series of Spherical Harmonics. 



The results obtained in this investigation lead to a new method of calculating the 
coefficients of the series of spherical harmonics which represent a function F of 
spherical co-ordinates. If the values of this function are given for all points of a 
sphere, the coefficients of the term involving T." is known to depend solely on the 
integral j'FT^ d$ taken over the surface of the sphere. This theoretically perfect 
proceeding has to be modified in practice when the values of F are known only at 
definite points, so that the intermediate positions have to be evaluated by a process 
of interpolation on the supposition that F is continuous everywhere. 

All methods which have been applied so far, suffer from, the serious inconvenience 
that the method of least squares is applied in such a way as to make the value of 
the coefficients of lower degrees depend on the number of terms which are taken 
into account. That is to say, the coefficients are not independent of each other as 
they ought to be. 

A. method suggested by F. E. Neumann, which is free from this defect, introduces 
other complications, and has never, as far as I know, been applied in practice. 

The theorem of § 8 offers a simple solution of the practical difficulties, and reduces 
the whole problem to an expansion by means of Fotjriek's analysis, which can be 
carried out either by the well-known process of calculation, or by mechanical means. 

Let F be expressed, in the first place, for different circles of latitude as a series 
proceeding by cosines and sines of the longitude and of its multiples. This first 
step is common to all methods. 

The result may be expressed symbolically by 

F = k {) + K f cos r/> -f- k~ cos 2(f) + . , , + K/ sin <f> -f K 2 sin 2c/> -f . . . , 

where k°, k\ K 7 , K. 3 . . . are functions of the colatitude. If their values are known 
at q ~* 1 equidistant circles of latitude, and at the poles, we may determine the 
coefficients a* and a", which satisfy the equations 

k* =s ciq + a? cos + a? cos 20 + . . . af, cos p0 + . . . a* cos q0, 

K " = < + a? cos 6 + a? cos 20 + . . . a* cos p0 + . . . a£ cos qd. 

If we give to the successive values 0, rr/q, 2rr/q . . . qrr/q, we shall have q + 1 
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equations to determine q + I coefficients. The solution is conducted according to 
well-known rules, making use of the proposition that, if p and s are integer numbers 
smaller than q, 



— i + cos -7 cos -— + cos ■■- - cos +*« 






, (? ™~ 1)^7T (# — 1)S7T . . 

+ cos 1 --~ cos + -A- cos prr cos am 

If jd and 5 are equal to each other the sum on the right-hand side is equal to \q, 
The factor ^ of the first and last term should be noted. If we designate by 

k , fc 1? k 2 . . . k p the value of k for # = 0, and at the successive circles of latitude, it 

follows that 

~~- aj = I (/Co + K 7 cos JP^) + 2 < cos spw/q. 

The condition under which the coefficients are determined is that the coefficients of 
the Fourier series higher than a\ are zero. It will be noted that, even if we suppose 
some of the lower coefficients to vanish, the equations will still give those values for 
the remaining coefficients which, according to the method of least squares, fit in best 
with the assumed values of k, but that in the calculations half weight only is given 
to the values at the two poles. 

[The above method of obtaining the coefficients a£, which is identical with that 
in common use, when the range of 9 is 2tt, is convenient whenever the function to be 
analysed is known at every point of the sphere, so that the coefficients k and K may 
be determined for a sufficient number of equidistant circles of latitude. But other 
methods are available, and hence the process of obtaining the coefficients of the 
series of spherical harmonics which I am endeavouring to explain, is not restricted to 
cases where the original function is known everywhere ; provided it is continuous, 
as well as its derivatives over the surface of the sphere. Graphical interpolation 
may be employed to determine the function at unknown points with sufficient 
accuracy, and there are several good mechanical devices in existence, by means of 
which at any rate the lowest and most important coefficients of the Fourier series 
may be found. It will, in some cases, materially help this process of interpolation if 
it is remembered that continuity at the pole involves the vanishing of all the values 
of k and K except k .— August 2, 1902. 

Having calculated the coefficients, the reduction to spherical harmonics is made, by 
substitution of 

cos pd = A^Q* + A.U1QU1 + ■ - • M-iQ^-i + • - • A^Q: + . . - , 

where the values of A^ may be calculated and tabulated once for all. By the 
theorem of § 8, all coefficients vanish up to and inclusive of A£_ 2 if a be odd, so that 
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the spherical harmonic of degree n will only depend on the Fourier coefficients 
Op, for which p is equal to or smaller than n + 1. When cr is even we secure the 
same advantage by developing k* and K°" in terms of the sine functions, and we 
write in that case : 



K 



°" = hi sin 6 + &? sin 20 + . . . &£ sinpO + • • • 



K°" = /^sin + yS^sin 2# + . . . ft* sin p/? + . . . 
and 

sin p0 = B^Q^ + Bj +1 Q; +1 + . . . B;Q; + . . . 

To calculate finally a coefficient such as B^ we proceed in the usual way, thus 



W 



Uherefore 



c 



^ 2??, + 1 (™ — O" — 1) M (^ — O) ! ! -.yr- , . .. 

B " = - — - ■--/• v ----; ttIVu when n is even and cr m even 

2 O + o-) ! ! (* + 2) ! ! 

= - — - - ■-<-■■■->- -- '- --- 7tOh when ?i- is odd and a* is even. 

2 (n + o- - 1) ! ! (n + 3) ! ! 

Similarly 

4 2'// + 1 (n — cr) l\ (n —■ o) ! ! TT , . . 

A°" = — - -/_ -A— — ■-—--, ttV^ when w, is even and cr is odd 

2 (w, -f <r - 1)!! O + 2)!! 

2% -f 1 (^ — o" ~~ -1) • i (^ — 2) ! ! TTn . -, * 1 i t*ii 

= ■- v — /__>_ _-/__ Tj-Lf ^ when n is odd and cr is odd. 

2 (ti + o-) ! ! (7i + 1) ! ! 

Combining equations, we find that if gl denote the coefficient of l^Ql cos crc/> in the 
levelopment of F, and -h? the coefficient of t%Q,l sin crc/>. 



? 



cf = ■— — /_>.- £- — 77- X (w} a£ V J ( r>), when n is even and cr odd and p odd 

= — — :: - -- ---- ^-i^_l_— A*: ^ ^(2) c^U^ ( p\ when ft is odd and cr odd and 73 even 

— — v y — v - J — it l*^ uZM." (p) when n is even and cr even and p odd 

2*S (76 + cr) 11(91+ 2)!! p = 1 p " u; l 

2?l + 1 (^ - CO ! ! (??, — 4) ! ! 2 Zm l 7 rraa- / \ t. 11 1 J 

= — — — ■ '--—-7 ~tt"~ ; ?r z ( } bZS* ( 7)) when n is odd and cr even and p even. 

2^ + <r — 1)!! (u + 3)!! v =9 2 xt } r 

To obtain h°„ substitute a and ft for a and b. 

The coefficient ^ is introduced because the quantity here designated by Q^ is not 
uniformly accepted as the standard form for a spherical harmonic. For some purposes 
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% 



it is more convenient to take the form of Gauss, for which Thomson and Tait have 
adopted the symbol 0£, and which is connected with Q* by the relation 



{Hj 



(ll •— or) : 






1 (2ra - 1) ! ! % " 



Another constant, which has been applied by Adolf Schmidt, is based on the 
consideration that for numerical work it is inconvenient to deal with functions the 
average values of which differ considerably from each other. This author therefore 
takes a function rj@; as basis of calculation, and determines r% so that the average 
value of the square of rl®% cos <x^> over a sphere of unit radius is equal to unity. The 
function B,£ defined in this way is connected with B" h and Q* by 

R' = (2n -IV \/-~-- A2n + l) 9"- = A/9J 2n ± 1 U n ~«) [ (> 
ix, n _ ^n, i).. v (% + ff) , (M _ ff) , u B - y (n + a) . V,. 

where e - is equal to 1 for cr = and equal to 2 for all other values of cr. 
If a function is to he expanded in terms of It;", we must therefore write 



or 



e^ (2n + 1) (?& — or) I 
(n + cr) ! 



For numerical work the introduction of li* in place of Q* j)ossesses undoubted 
advantages. Although at first I was reluctant to adopt the additional complication 
due to the introduction of a square root and the addition of yet another function to 
those given by previous writers, I found that the inconvenience of tabulating values 
differing considerably in magnitude from each other was very great, and I therefore 
felt myself almost compelled during the course of the investigation to adopt Schmidts 
function as above defined. 

We now substitute t" t into the equations for (f lh and h^ and write : 



m (n -3)1! /'2n + 1 fi / ln~+~~crj\\(ii-cr))l 



a 



X7V 



Vn 2\*>~ +'->)]]' V €(f 'V ( ?l + £r _l)!J( u _ <r -.l)!! ' 



7T (n - 2)11 /'In + 1 7(72, + a- - 1) !J_(?t - cr - 1) ! ! 
Y (??•"+" 1)YI "V 6<r " V "( n + -)!!(n - cr) ! f """ 



7T 0/, - 3) ! ! /2?fc + 1 /O + cr - 1) ! ! (w - a- - 1) ! 



™* — - >— .~ --'-.' a / —:....: : a / \----. v .....-' ^ _ J _ ivi ff 



7T (-;?. - 4) ! ! /2u + 1 / (■«, -f- a-) ! ! (■>*, — cr N f ' 



71 ~~ 2 (?* + 3) ! ! V 6<r V ( 7l + o- _ 1) ! ! ( ?l _ ff _ 1) ! ! °»' 

where c = 1 and e^ = 2 if cr > 0, 

By substitution of the values of V£, &c,, we may also write more simply 
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+ 1 



L 'i I 



^£<r ► -1+1 



Hi cos pOdfi ; if a be odd, p odd, n even. 



ul = -— I II* cos p0dfi ; if cr „ „ p even, n odd. 



9*? i 



m cr _ i j^o- g j n p^ ^ . j£ ^ | )e even? p 0( j ( | ? n even. 



•"-(T •> — 1 



I r + i 

5 (r = __ p^ smpddfi ; if or „ „ p even, n odd. 

It would thus appear to be unnecessary to have separate symbols for u" and if a or 
for ml and si, but as it is convenient to tabulate separately the integrals for odd and 
even values of p y the retention of all four symbols facilitates reference. 

We may now obtain the final coefficients by summation thus : 

p — n + 1 

gl = %® a^;^, when n is even, cr odd, and p odd, 

2)=:'/) + l 

= S (2) a^4, „ n„ odd, cr „ „ p even, 

2> = '/l + l 

= S (2) &£ra£ ? „ n ,, even, cr even, „ p odd, 

^— -,? + i 

= £ (2) &£<, ?, w „ odd, a* „ ,, p even, 



0) — *> 



with similar equations for /*£, a and h being replaced by a and /3. - 

The values of vl, ul, ml and si are given in Tables V.— VIII. and their logarithms 
in Tables IX. — XII. , for values of n, cr and p up to 12 inclusive. By means of these 
tables we may, for instance, write down at once the coefficients as far as the third 
degree as follows : 

q° = -6801756o° 

#/ = - •2195256/4- -6585756./ 

g° — - '25974G6/ 3 + 'G49365&.,, 

r// = -6801.750/ - '340087a/ 

<■// = -380229a/ — -380229a,' 

gj = -159061a/ + -318123%' — -397653a./ 

a,/ = -5703456/ — -1901156, 2 

gf — -5029966/ - '2514986/ 

g/' = -616042a,, 3 '— '410695a, 3 -f •102672a./. 
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The calculation of the ten coefficients g up to the third degree involves therefore the 
determination of sixteen Fourier coefficients. There is no h coefficient of zero type, 
and there are therefore only six lb coefficients up to the third degree which depend 
upon twelve Fourier coefficients. In general there will be -| (n '+ 1) (n + 2) 
g coefficients, and -J(n + 1) ft? h coefficients up to and including the n th degree, giving 
a total of (n + I) 3 coefficients in the series of spherical harmonics. To determine 
these completely;, it would be necessary to calculate (n + l) 3 Fourier coefficients 
for the quantities g and n n + 1 for the quantities h, giving a total of (n + 1) ( n + 2) 
Fourier coefficients. The tables accompanying this paper being calculated for 
values of n up to 12 are therefore sufficient to determine 169 coefficients in the 
series of spherical harmonics, by means of a simple computation after the deter- 
mination of 182 coefficients of the Fourier series. These latter coefficients may be 
evaluated by mechanical devices, 



§ 12. Special Application to the Theory of Terrestrial Magnetism. 

The advantages of the proposed method of obtaining the coefficients of a series 
of spherical harmonics are considerably increased, when the quantities to be repre- 
sented by such a series are not given directly, but by means of their differential 
coefficients. This is the case when the magnetic potential has to be calculated by 
means of the observed magnetic forces. If X and Y represent the components of 
magnetic force resolved towards the geographical north and east respectively, the 
magnetic potential is determined by means of the equations 

dV/dd = X ; dV/d<f> = - Y sin 0. 

If no electric currents of sufficient intensity traverse the earth's surface, a 
function Y can be found which satisfies both equations. If Y sin be obtained in a 
series proceeding by spherical harmonics, then all the terms which depend on the 
longitude <f> are at once expressed in a similar series, as the integration according 
to <f> leaves each term in the standard form. The treatment of the other component 
acting along the meridian involves, however, serious difficulties, and it is not 
necessary here to enter into the question as to the more or less complicated methods 
by means of which Y has been hitherto derived in the standard form from X. 

The method based on the results of the preceding investigation avoids these 
difficulties. For odd values of <x, X is expressed in a series, each term of which 
has the form cos n0 sin m<$> or cos n0 cos m<£ ? while for even values of cr, X is 
expressed in terms of the form sin n0 sin mcf) and sin n0 cos m<£. After integration 
with respect to 0, the formula of § 11 will determine the required coefficients. We 
may treat the eastern force similarly, obtaining a series proceeding by cos^;>0 or 
sin p0 according as cr is odd or even. Y sin is then derived in a series proceeding 
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by sin pd or cos pd, and after integration with respect to (j> the equations of § II 
may be applied. The spheroidal shape of the earth may also, if necessary, be easily 
taken into account, but I reserve the discussion of this matter until the completion 
of some calculations, on which I am at present engaged, will allow me to compare 
the results obtained by the method I advocate with those obtained in other ways. 

Great pains have been taken to guard against numerical errors, and it is hoped 
that all numbers given in the tables are correct, with the proviso, however, that 
the last decimal place in Tables V. and XII. is uncertain by one unit, and wrong 
possibly even occasionally by two units. Professor Core has assisted me in the 
numerical work. 
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Table L— Values of V~n(p)- 



cr=l 






o-^O 



fT 



7 



cr -. 



(T - 



= 11 



w.; p\ 



2 



6 

8 

10 

1° 

J. w 



4 
6 

8 
10 
12 



6 

8 
10 
12 



8 
10 

1 o 

1 *J 



10 






J. i«J 



1. 



1 
1 
1 
1 
1 
1 



3 
3 



3 



3 



o 







5 



7 
7 
7 



9 



9 



11 



3. 



1 



7 

5 

6 



4 y 

1 ft 

l 2' 



9 

~2 

1 



9 



117 

TO" 

ft 1 

Til" 



-9 



2 



4 
.1- <y 



-H 



4 1 

7 



i ft r> 



4 i 

' 4"" 



9 u 

' 4 



•). 



2 1_ 

5 

« 

TV " -*" 
■J t> 

1 4 ft 



.« %' ^ 



3 

98 

1 7 9 
10 8 



5 
10 

8 ft 

1 4 






10 



1-1 5 

"IT" 

1 ii 1 

1^ 



' 1 1 " 



ii t> 



7. 



.. 3 .*L 

420 

To" 
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